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Abstract 

In this paper we will discuss some features of the bi-Hamiltonian 
method for solving the Hamilton-Jacobi (H-J) equations by Sepa- 
ration of Variables, and make contact with the theory of Algebraic 
Complete Integrability and, specifically, with the Veselov-Novikov no- 
tion of algebro-geometric (AG) Poisson brackets. The bi-Hamiltonian 
method for separating the Hamilton-Jacobi equations is based on the 
notion of pencil of Poisson brackets and on the Gel'fand-Zakharevich 
(GZ) approach to integrable systems. We will herewith show how, 
quite naturally, GZ systems may give rise to AG Poisson brackets, 
together with specific recipes to solve the H-J equations. We will then 
show how this setting works by framing results by Veselov and Pen- 
sko'i about the algebraic integrability of the Volterra lattice within the 
bi-Hamiltonian setting for Separation of Variables 



1 Introduction 

The Hamilton-Jacobi equations, and the problem of their separability, are one 
of the many fields in mathematics in which the influence and heritage of Carl 
Gustav Jacob Jacobi is still alive. Such a problem, which can be considered 
one of the fundamental problems of Theoretical Mechanics, is rooted in the 
foundational works of Jacobi, Stackel, Levi-Civita and others. It has recently 



1 



received a strong renewed interest thanks to its applications to the theory of 
integrable PDEs of KdV type (namely, the theory of finite-gap integration) 
and to the theory of quantum integrable systems (see, e.g., fo\ l2*Uj). 

The constructive definition of separability originally due to Jacobi is the 
following. 

Let us consider an integrable Hamiltonian H on a 2n-dimensional phase 
space, that is, let us suppose that, along with H = H\ we have further n—1 
mutually commuting integrals of the motion H 2 , . . . H n , with dHi A ... A 
dH n ^0. 

Definition 1.1 An integrable system (Hi, . . . , H n ) is separable in the canon- 
ical coordinates (p, q) if there exist n independent relations 

$ i (q i ,p i ;H 1 ,...,H n ) = 0, i = l,...n, (1.1) 

connecting single pairs (qi,Pi) of coordinates with the n Hamiltonians Hj. 

The link of this definition with the theory of those integrable systems that 
admit a Lax representation with spectral parameter 

L(\) = [L(X),P(X)), 

such as those associated with classical limits of quantum spin systems and/or 
those coming from suitable reductions/restrictions of KdV-like evolutionary 
PDEs, is self evident. Indeed, the Lax representation of a system provides 
us with a natural candidate for the separation relations: the characteristic 
polynomial of L(\), also known as the spectral curve associated with L(X). 
However, the possibility of successfully applying the Lax method relies on 
three non-algorithmic steps to be solved: 

• To find the Lax representation of a dynamical system; 

• To prove that the spectral invariants of L(X) are mutually in involution 
(i.e., to prove that the Lax representation is compatible with a classical 
r-matrix structure); 

• To give canonical coordinates as the coordinates of n points lying on 
the spectral curve, i.e., to actually implement what is sometimes called 
Sklyanin's magic recipe. 

The setting devised by Veselov and Novikov [22| to characterize algebraic 
integrability requires that the phase space A4 of a Hamiltonian system fulfill 
some properties. They can be, quite roughly, summarized as follows: 
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a) M. has the fibered structure 



M 



B 



(1.2) 



where the base B is a n-dimensional manifold whose points determine 
an algebraic curve T(b), and the fiber is the k-th symmetric product 
of that curve. In more details, one requires that T(b) be given as an 
m-sheeted covering r(6)-^C of the complex A-plane, and that points 
of M. can be parameterized via the curve T(b), and a set of k points on 
it, that is, the coordinates Ai,. . . ,\k of the projection on the A-plane of 
a set of points on it, as well as discrete parameters q that specify on 
which sheet of the covering the points live. 

b) An Abelian differential Q(T) on T (or possibly on a covering of F), 
smoothly depending on the points b G B, is defined. It is furthermore 
required that, if Q(T) is given by 



according to the representation of T as a covering of the A-plane, the 
closed two-form 



give rise to a Poisson bracket, conveniently called algebro- geometric 
Poisson bracket, with A* and /x 4 = Q(b; A*) playing the role of Darboux 
coordinates on the symplectic leaves of this bracket. 

In such a case, it was proven in |22] that functions that depend only on 
the curve T - i.e., on the points of B - are in involution with respect to the 
Poisson bracket defined by (|1.4|) . and these geometric data explicitly define 
action-angle variables for these Hamiltonian flows. In that fundamental pa- 
per it has also been shown that a number of integrable systems, of classical 
(i.e., mechanical) type as well as obtained by suitable reductions of soliton 
equations, can be framed within such a scheme. In particular, the paper [23] 
shows how the Volterra lattice fits in it. Finally, it is worth mentioning that 
Sklyanin's method [20] of the poles of the Baker- Akhiezer function, originally 
introduced in the study of Hamiltonian systems as a byproduct of "Quantum 
Integrability" , can be seen as a particularly efficient scheme of implementing 
the Veselov-Novikov axiomatic picture. 



Q(T) = Q(b; X)d\ 




k 
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More recently, a bi-Hamiltonian approach to Separation of Variables 
(SoV), has been exposed in the literature (see, e.g., [SI El El UH1 EEU 12] )• 
Such a scheme can be seen as a kind of bridge between the classical and the 
modern points of view, putting an emphasis on the geometrical structures 
of the Hamiltonian theory. In this framework it is possible to formulate in- 
trinsic conditions on the integrable system (Hi, . . . , H n ) to a priori ensure 
separability in a set of canonical coordinates. It requires the existence, on 
the phase space (M,u>), of a second Hamiltonian structure, compatible with 
the one defined by cu. Namely, the bi-Hamiltonian structure on M allows, as 
it has been shown in a number of examples, 

1. To encompass the definition of a special set of coordinates, to be called 
Darboux-Nijenhuis (DN) coordinates, within a well defined geometrical 
object. 

2. To formulate intrinsic (i.e., tensorial) conditions for the separability 
of a Hamiltonian integrable system, in the DN coordinates associated 
with the bi-Hamiltonian structure. 

3. To give recipes to characterize, find and handle sets of DN coordinates. 

In particular, in [TT. a detailed discussion of the bi-Hamiltonian scheme for 
SoV in the case of Gel'fand-Zakharevich type jTH] systems was presented. It 
was also pointed out how the separation relations of such systems were, under 
genericity assumptions, of degenerate type, that is, the functional form of the 
separation relations (jl.lj) is the same for all pairs of separation coordinates 
(Aj,/ij), which essentially means that these coordinates are the coordinates 
of different points on the same algebraic curve. 

In this paper we want to elaborate further on this issue, and, in particu- 
lar, establish a connection between the bi-Hamiltonian scheme and the VN 
setting. This will be done in the first part of the paper, and, namely, in 
Proposition 13.31 which shows that (under suitable assumptions) DN coordi- 
nates associated with Gel'fand-Zakharevich systems can be seen as algebro- 
geometric canonical coordinates in the VN sense. In the second part of the 
paper we will apply our scheme in revisiting the algebro-geometrical integra- 
bility ^3 EI] of the well-known Volterra lattice. 

More in details, the paper is organized as follows. In Section |21 we will 
briefly review the bi-Hamiltonian set-up for SoV. In Section we will recall 
some results exposed in ^T] , and discuss the relations of the bi-Hamiltonian 
approach with the VN scheme. In Section HI we will collect a few results 
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concerning the algebro-geometric scheme of integrating the Volterra lattice. 
Finally, in Section |S] we will show how the bi-Hamiltonian picture of Section 
121 can be successfully applied to the lattice, with no significant differences 
between the cases with odd (resp., even) number of sites. 

2 Bi-Hamiltonian geometry and Separation 
of Variables 

The basic geometrical notion underlying the bi-Hamiltonian scheme for sep- 
aration variables is that of "semisimple ujN manifold". An uoN manifold is 
a symplectic manifold (M,u>) endowed with a second (possibly degenerate) 
Poisson tensor Pi which is compatible with the Poisson tensor P associated 
with the symplectic form u. This means that P\ = Pi — XPq is a Poisson ten- 
sor for all AeR. In this case P\ is called the Poisson pencil and (M, P ,Pi) 
is a bi-Hamiltonian manifold. It can be shown (see, e.g., [T5] ) that the (1,1) 
tensor field N = Pi o Pq 1 has the property 



for all vector fields X, Y on M, that is, its Nijenhuis torsion vanishes. The 
tensor field N is called Nijenhuis tensor or recursion operator of the uN 
manifold. It turns out that the characteristic polynomial of iV is the square 
of a polynomial A (A); the uN manifold M is called semisimple provided that 
the roots of A(A) be (generically) simple. 

A special class of coordinates, to be called Nijenhuis coordinates, are 
provided by the spectral analysis of the adjoint recursion operator iV* = 
Pq 1 o Pi. Indeed, one has the following results (see, e.g., |TH llb[ Ho]): 

1. The eigenspace Aj corresponding to any root A, is an integrable two- 
dimensional codistribution, that is, one can find n pairs of functions 
fi,gi (to be called Nijenhuis coordinates) satisfying 



2. The eigenspaces Aj and Aj are orthogonal with respect to the Poisson 
brackets induced both by P and Pi. This crucial property can be 
very simply proven. Indeed, let / and g be such that their differentials 



[NX, NY] = N([NX, Y] + [X, NY] - N([X, Y})) 



(2.1) 



N*dfi = Xidfi, N*dgi = \dgi\ 



(2.2) 
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belong respectively to Aj and to Aj, with i ^ j. Then one has, on the 
one hand: 

{f,g} x = (df,P l( lg) = (N*df,P dg) = \{f,g} Q . 

Switching the role of / and g one sees that {/, g}\ = \j{f, g}o, whence 
the assertion. 

3. Since (as it is easy to prove) the Poisson bracket {fi,gi}o with respect 
to Pq (and also the one with respect to Pi) of functions that satisfy ()2.2|) 
still satisfies N*d{fi, gi} = \{fi, gi}o, it is possible to parameterize Aj 
with a set of coordinates Xi,yi, called Darboux-Nijenhuis coordinates, 
that are Nijenhuis coordinates and are canonical for uj (whence the 
addition of "Darboux" in their denomination). 

4. It is important to notice that the Nijenhuis tensor A" of an ujN manifold 
defines, at each point m G M, a linear operator N m : T m M — > T m M. 
As such, its eigenvalues (that are point- wise the roots of det(A m — A)) 
may depend on the point m. In this case we will call these roots, 
nonconstant roots of A(A). If A is a nonconstant root of A(A), then it 
satisfies the characteristic equation 1 ()2.2j) . 

N*d\ = XdX 

The following proposition has been proven in [TT] . 

Proposition 2.1 Suppose that (M,u,Pi) is a semisimple uN manifold of 
dimension 2n. Let Xi, . . . , x n , yi, . . . x n be Darboux-Nijenhuis coordinates on 
M and let F 1; . . . , F n be functionally independent Hamiltonians, that are in 
involution with respect to the Poisson brackets induced by Pq and P\. Then 
the Hamilton- J acobi equations associated with any of the Hamiltonians Fi 
can be solved by additive separation of variables in the Darboux-Nijenhuis 
coordinates {xi,yi)i 

=l,...,n ■ 

The foliation given by the functions Fj will be called a bi-Lagrangian folia- 
tion. Such foliations provide a geometrical description of separable systems, 
exactly like Lagrangian foliations describe integrable systems. 

To elaborate further on the geometric structure of uN manifold, it is 
convenient to suppose that the eigenvalues of the Nijenhuis tensor A^ be 

Actually, the constant eigenvalues trivially satisfy the same equation. 
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functionally independent. A straightforward observation [Tl^ ITT] shows that 
one can compactly write the characteristic equation N*dX{ = XidXi in terms 
of the minimal polynomial 

idim(M) 

1=1 

as the polynomial relation 

N*dA N (X) = XdA N (X). (2.3) 

Actually, relations of this kind are very important for our purposes. Indeed, 
in JTJ we proved the following proposition: 

Proposition 2.2 Let $(A) be a smooth function defined on theuoN manifold 
M, depending on an additional parameter X. Suppose that there exists a one- 
form a$ such that 

N*d$(X) = Ad$(A) + A N (X)a^ . (2.4) 
where A^ is the minimal polynomial of N. Then: 

a) the n functions $j obtained evaluating the "generating" function $(A) for 
X = Xi,i = 1, . . . ,n are Nijenhuis functions, that is, they satisfy N*d&i = 
Xid$i. 

b) J/$(A) satisfies \2.Ji\j and Y\ = —P^dpi are the vector fields associated via 
Pq to the coefficients of the minimal polynomial of the Nijenhuis tensor, then 
all functions 

^(X) = LieyXHX))) 

satisfy \2.4\j as well. 

c) In particular, if $(A) satisfies, along with \2.$ , the relation 

Lie Yl (<$>(X)) = 1 modA 7V (A), (2.5) 

then the functions (Aj,/x, := 3>(Aj)) provide a set of Darboux- Nijenhuis coor- 
dinates on M. 

Definition 2.3 We will call a generating function $(A) satisfying equa- 
tion l[2.4\ ) a Nijenhuis functions generator. 
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3 Gel'fand-Zakharevich systems 



In many of the models considered in the so-called "Modern Theory of In- 
tegrable Systems", that is, finite-dimensional Hamiltonian systems obtained 
as reductions of integrable PDEs, and/or classical analogs of quantum spin 
systems, two instances occur: 

• The phase space M of the system is not the cotangent bundle to a 
smooth manifold Q. 

• M is endowed with a pair of compatible Poisson brackets, but none of 
them is nondegenerate. 

This is, for instance, the case of the Volterra lattice. This geometrical in- 
stance has been formalized in a series of papers by Gel'fand and Zakharevich. 
In particular, under some technical assumptions, one of which is that the 
corank of P equals that of Pi, it is possible to construct iV = corank(P ) 
Lenard-Magri sequences that start with a Casimir function of Pq and end 
with a Casimir function of P\. The Hamiltonians of these sequences can be 
conveniently collected in polynomials H^ a '(X), a = 1, . . . , N, in the variable 
A, satisfying 



called polynomial Casimirs of the pencil. It is nowadays customary to denote 
such a geometrical instance of bi-Hamiltonian manifolds (M, P , Pi) with the 
name of Gel'fand-Zakharevich (GZ) manifold. 

The collection of the degrees of the polynomial Casimirs is a numeric 
invariant of the bi-Hamiltonian manifold. Notice in particular that Casimirs 
of degree are nothing but common Casimir functions of the two Poisson 
tensors. We will, in the sequel, call these common Casimirs trivial Casimirs, 
and refer to the others (namely, those originating a non void Lenard-Magri 
sequence) as nontrivial ones. 

In the case of bi-Hamiltonian systems associated with (coefficients of) 
polynomial Casimirs, one can try to use the bi-Hamiltonian scheme for SoV 
reducing the systems to a suitable uN manifold. More precisely, one can 
consider a symplectic leaf S C M of Pq and a suitable deformation of the 
Poisson pencil, discussed in detail in (TU1 HH [T7j . 

First one fixes a maximal set Ci , . . . , of independent nontrivial Casimirs 
of Pq, and finds k (independent) vector fields Z\, . . . ,Zk such that 




(3.1) 



Lie Za (C b ) = 5 ab 



Lie Za K a = 
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where K a are the common Casimirs 2 . We will hereinafter refer to the distri- 
bution Z generated by the k vector fields Z^ as the transversal distribution. 

Then one considers the vector fields X a = PidC ai for a = 1, . . . , k, and 
the "deformed" tensor 



that restricts to S. If the algebra of functions vanishing along Z%, . . . , Z^ is a 
Poisson algebra for the pencil Pi — APo, it turns out jTUj that the deformation 
Pi defines on S a Poisson tensor P\\S compatible with the restriction Po\S 
of P to its symplectic leaves. So S is endowed with the structure of a uN 
manifold. By the definition of Pi, the restrictions of the coefficients 
of P^ a )(A) to S will be separable Hamiltonians in the Darboux-Nijenhuis 
coordinates defined on S (provided S is semisimple). 

Remark. To summarize, the ideas underlying such a reduction pro- 
cedure are the following: the Gel'fand-Zakharevich scheme provides - un- 
der some technical assumptions - a way for defining, via Magri-Lenard se- 
quences, a distinguished integrable distribution A on a bi-Hamiltonian man- 
ifold (M, P , Pi). It is called the axis of M, and is generated by the Hamil- 
tonian vector fields associated with the coefficients of the Casimir poly- 
nomials, so that the leaves T of A are defined by the requirement that these 
coefficients be constant along A. 

If we fix our attention on one of the elements of the Poisson pencil, say, 
Po, and consider its symplectic foliation S, we have that T H S defines, on 
the generic symplectic leaf S of <S, a Liouville integrable system. 

However, since in general the symplectic foliation associated with the 
other Poisson tensor Pi does not coincide with S, S does not come equipped 
with a natural bi-Hamiltonian structure. Finding the distribution Z with the 
properties outlined above amounts to finding a deformation Pi of Pi that 

a) endows the symplectic manifold S with a compatible second Poisson 
tensor, and hence with the structure of a uN manifold. 

b) Preserves the commutativity of the Hamiltonians H\ ^ , that is, provides 
T fl S with the structure of a bi-Lagrangian foliation. 

In the papers referred to above, the distribution Z spanned by Zi, . . . , Zk was required 
to be satisfy the stronger condition T p M — T p S © Z p for all p € S. Actually, as it should 
be clear form common Casimirs do not enter the reduction procedure. A similar 
instance with common Casimirs has been considered in j!2j . 
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(3.2) 



a=l 
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Although to find the transversal distribution Z is a non algorithmic pro- 
cedure, we notice that this is a quite efficient way of providing the symplec- 
tic leaves S of Pq with a compatible Nijenhuis tensor (a problem which, in 
principle, requires the solution of a system of nonlinear partial differential 
equations). 

Definition 3.1 We say that a GZ manifold (M, Pq, Pi), endowed with a 
transversal distribution Z satisfying the above mentioned assumptions, ad- 
mits an ajfine structure if it is possible to choose a complete set of nontriv- 
ial Casimirs of Pq, and a corresponding basis of normalized flat generators 
{Zb} b=i,...,k i n Z such that, for every Casimir of the Poisson pencil H a (X) 
and every b, c, one has the vanishing of the second Lie derivative of the 
Casimir polynomials: 

Lie Zb Lie Zc (H a (\))=0. (3.3) 

The above definition might seem somewhat ad hoc. Its relevance can be 
summarized in the following points (whose proof can be, once more, found 
in jTTj) that hold in the case of affine structures. 

• The nonconstant roots Aj of the determinant of the matrix whose entries 
are 

g ab = Ue Zb {H a {\)) 

satisfy 

N*d\i = XidXi, 

namely, they are roots of the minimal polynomial of the Nijenhuis ten- 
sor induced on (any of) the symplectic leaves of P by the pair (P , Px). 
In particular, if there is only one non-trivial Casimir polynomial H(\), 
and Z is a corresponding normalized flat generator, the nonconstant 
roots of the polynomial Lie^(if(A)) are "nonconstant" eigenvalues of 
the Nijenhuis tensor N. 

• The separation relation satisfied by the non-trivial Hamiltonian func- 
tions and the Darboux-Nijenhuis coordinates are linear in the Hamil- 
tonians, that is, are of (generalized) Stackel type. 

The following Lemma, whose proof is a simple application of some notions 
of Poisson geometry, will be frequently used in the sequel. It provides a link 
between the properties of functions on the GZ manifold M, which depend 
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polynomially on the parameter A of the Poisson pencil P\ — XPq defined on 
M, and the properties of the evaluation of such functions in A = \ w.r.t. 
the induced Nijenhuis structure on the symplectic leaves. In plain words, it 
allows us to work on the GZ manifold M, without having to actually perform 
the reduction procedure. 

We still suppose that (M, Po, Pi) is a GZ manifold, with k non-trivial 
Lenard-Magri sequences. We suppose that Zi, . . . ,Zf~ are normalized transver- 
sal generators for the distribution Z we considered above. We recall that, in 
this situation, the symplectic leaves of P are uN manifolds, with induced 
Nijenhuis tensor N = Pq P\.. 

Lemma 3.2 Let F\ be a function on M , invariant along the fields Z i; that 
defends holomorphically (say, polynomially) on an additional parameter A; 
its restriction f x to a symplectic leaf S of Po satisfies the "eigenvector" equa- 
tion 

N*df x \ x=Xi =\ t dfx\ x=Xt 

for all eigenvalues A, if and only if the following equality holds, parametrically 
in X, on the GZ manifold M: 

k 

{G,F x } Pl -J2^e Za (G){K[ a \F x } Po = \{G,F} Po , (3.4) 

a=l 

for any G G C°°(M), where the K[ a) satisfy P x dC a = P dK[ a \ a = l,...,k. 
Otherwise stated, we have to require that 

T\dF x = XP dF x . 

where Pi is defined in \3. fy) 

A direct generalization of the above proposition shows that spectral curves 
might be a source for finding Nijenhuis functions generators. 

Proposition 3.3 Let us consider a generating function r(A, fj,) of Casimirs 
of a Poisson pencil P x , and let us suppose that T(A, /i) — defines a smooth 
algebraic curve. Let S be a generic symplectic leaf of Po and let N be the 
Nijenhuis tensor associated — according to the scheme outlined above — with 
P x and a suitable transversal distribution Z. Suppose that f is a Z -invariant 
root of the minimal polynomial of N, i.e., 

N*df = fdf, and Zi(f) = 0, i = 1, . . . , k, (3.5) 
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and suppose that r(/x, /) = defines generic point(s) of the affine curve 
T(A,/i) = 0. Then, any solution g of the equation T(g,f) = which is 
invariant as well under Z satisfies N*dg = fdg. 

Proof: We first notice the following. Let us consider a bivariate polynomial 
F(\,[jl) = Ylij f(ij)^ 1 ^ i with coefficients f^j) that are functions defined on 
a manifold M, and two more distinguished functions on M, say / and g. If 
we define T := F(f,g), then: 



dT = dF(X, n) 



dF 

M=g OX 



df + — 

M=9 dfl 



-.fdg. 



(3.6) 



We consider the equation of the "spectral curve", r(/x, A) = 0; so we get for 
the zeroes of the function 1Z = T(f,g), 



= dT(\,fi) 



or 



^jdf + ^ 



M=9 



dg. 



(3.7) 



Let us suppose, for simplicity, that T(A, /x) = T (X{i) + Yli=i /U n 'iP(A), where 
r (/U, A) is a constant polynomial (possibly depending on the common Casimirs), 
and Ui ares suitable integers. Considering the action of the k transversal vec- 
tor fields Zi, we get the k equations: 



Lie^(^) 



dr 



dT 

, 4 + Lie Zi (/)- 



A»=9 



5> B 'Lie*(iP(A)) 



M=9 



0. (3.8) 



Applying — /P to eq. (|3.7|) we get, using Lemma I3~21 and taking into 
account that (Pi — XPo)dT(X, fi) = 0, 



8T 
dfi 



M=9 



i=l j=l 

™™ (3-9) 
Plugging into equations ()3.8I3.9|) the hypotheses on / we arrive at the system 



fr(/,s) = o 

Lie Zi ( 5 )— ^-^^Lie Zi (^(A)) 



A»=9 



0, i = 1, . . . , k. 



3=1 



(3.10) 



ar 
9/i 



A»=9 



(iV * _ j)^ _ ^^^Lie Zj /T(A)| A=/ d^ = 0. 

i=X j=l 
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The thesis follows noticing that if g is invariant under Z then this system 
reduces to 

r(/,<7) = o 
or 

and taking into account that, for T(A, //) smooth, the solution of system 



(3.11) 

^ f (N*-f)dg = 0, 

P-=9 



r(/,<7) = o 
or 



M=9 



<9/i 

are the (fixed) ramification points of T(A, //) = 0. 

□ 

This proposition provides us with the desired link between the bi-Ham- 
iltonian approach and the VN axiomatic picture of AG brackets. Indeed, it 
can be restated as follows: suppose we can find (by means of a Lax represen- 
tation, or by other means) a generating function for the Casimirs polynomials 
of an affine GZ pencil in the form of a bivariate polynomial T(X,/j,). Then, 
suppose that the coordinates of the points on the curve T(A,/x) = (whose 
A-projections give the roots of the minimal polynomial of the Nijenhuis ten- 
sor iV induced on symplectic leaves of Pq) satisfy the invariance condition 
specified in the above proposition. Then they are Nijenhuis coordinates, and 
so their Poisson brackets are given by 

{Aj,//j}o = Sij(pi(Xi, Hi), {Aj, fj,j}i = 5ij\i(pi(\i, Hi), 

Under the further assumption of irreducibility of the minimal polynomial of 
N, one sees that the unknown functions ipi cannot explicitly depend on the 
index i. Thus the formal integral 

will give the VN meromorphic differential defining Algebro-geometrical Pois- 
son brackets corresponding to P (as well as Q' = ^ gives those correspond- 

A 

ing to Pi). 
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4 The Volterra Lattice 

The Volterra lattice equations are the following set 

Q = - (4.1) 

which we consider to be defined on a periodic lattice q > 0, Q +n = q. They 
are generalization of the famous Volterra equations describing time evolution 
of competing species. 

The phase space of the VL can be seen as the restriction to the sub- 
manifold of vanishing momenta of the periodic n-site Toda Lattice. It is 
well-known HH1 123 ED Hj that eq. (|4.1|) are isospectral deformation equa- 
tions for the periodic difference second order operators of the form 

(£0)* = afc+i^fc+i + aki^k-i, a n+i = cti, ip n+i = Xipi. (4.2) 



where = yfck- In complete analogy with the Toda case, it admits a dual 
Lax representation 3 in terms of a 2 x 2 matrix C being given by the (ordered) 
product of site matrices 

£' = 44-1 W= ( \ q V 

The Volterra Lattice equations admit a bi-Hamiltonian formulation. In- 
deed, if one considers the quadratic Poisson tensor 

P ij = acjiSi+xj - 5 j+1>i ) (4.3) 

and the cubic one 

Q %3 = CiCj(ci + Cj)(5 i+ld - 8j +1>i ) + CiCi + iCi +2 5i + 2,j - c i Cj_ic i _ 2 5i-2,j (4.4) 
one notices that (|4.1|) can be written as 



where 



h 



j n n 

-logjjci, k = J2^. (4.6) 



t=i i=i 



3 This duality involves also, as in the case of the Toda Lattice, an exchange of the roles 
between the spectral parameter A and the eigenvalue fi. 



14 



We can collect those remarkable results by Penskoi and Veselov-Penskoi 19, 
123] . that will be used in the sequel, as follows. 

Let us set n = dllLi^) 1 > an< ^ consider the (normalized) characteristic 
equation 

T{X,n) = -Det(/i - £(A)) = 0. (4.7) 



7T 



Then it holds: 

1. T(A, n) = X + v — i^(/i), where the polynomial H(jj) is expressed as 



E- „(-i)V"- ° ' 7, 



2fc+l-2i^ 
7T 



/" • (ELoC-l)^ 2 ^^^), if n = 2A; + 1 



J; 



7T 



• (ELo^ 1 )*/^"^) + Wa+i, if n = 2A; + 2 

(4.8) 

The functions J, can be usefully characterized in the following way: Jj 
is the sum of all possible monomials ci x c\ 2 ■ ■ •q i of lenght z, where the 
indices l p are all different and not congruent to 1 modulo n. In |2"3~] 
subsets {Zx, • • • , k} C {1, . . .n} satisfying this property are called totally 
disconnected and we will adopt this definition in the next subsection. 
Notice that, e.g., J = 1, J% — Ei c «- Furthermore, notice that for 
n = 2k + 2 the last Hamiltonian T~tk+i is given by 



U 



fe+i 



A + — , with A 
A 



C2 C4 • • • C2fc+2 
Cl C4 • ■ ■ C2fc+1 



(4.9) 



2. (Theorem 3 of JH]). The functions Tii satisfy the Lenard-Magri recur- 
sion relations 

PdH = 0, PdUi = QdHi-ui = l,...,k, QdH k = 0. (4.10) 

For n = 2k + 2 the last Hamiltonian 7Yfc+i is in the kernel of both P 
and Q. 

3. Let {Ci, . . . , C/c} be a suitably chosen subset of k poles of a suitably nor- 
malized Baker- Akhiezer function \1/ associated with the Lax operator L 
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of (|4.2|) . and let Aj be corresponding coordinates on the spectral curve 
T. Then the coordinates {Q, Pi = 21o f } parametrize the symplectic 
leaves of P and satisfy 

{d, Pj}p = {&, Pj}q = (fhj. (4.11) 

We remark that one can compactly restate the results of item 2) in the 
following form: 

Proposition 4.1 Let us define the Poisson pencil 

P^ = Q- p 2 P, 

where P and Q are the quadratic and cubic Poisson tensors \4-%\4-4\) - Then 
the polynomial H(fi) defined by fl^.&| ) is a Casimir polynomial of the pencil 
PyP. . Moreover, for n = 2k + 2 the function 7ik+i is a common Casimir of 
the two basic elements P and Q of the pencil. Since functional independence 
of the Hamiltonians Hi is self-evident, we see that, in the GZ terminology, 
the phase space M of the n-site Volterra lattice, equipped with the Poisson 
structures P and Q, is a GZ bi-Hamiltonian manifold. If n — 2k + 1 or 
n = 2k + 2, we have a single non-trivial Lenard-Magri sequence comprising 
k vector fields. If n = 2k + 2 is even, the GZ manifold has a trivial Casimir 

5 Separation of Variables for the Volterra Lat- 
tice in the bi-Hamiltonian setting 

In this section we will make the final contact between the known results we 
collected above, and explicitly show how the picture of j2H] can be naturally 
framed within the bi-Hamiltonian theory of Separation of Variables for GZ 
systems described in Section |2] 

Our starting point are the GZ formulation of Penskoi's results, collected in 
Section 3, as well as the normalized spectral curve equation ()4.8|) . Obviously 
enough, we will set, as natural parameter of the Poisson pencil, the quantity 

v = [l 2 . 

We remark that both in the even and odd number of sites (or species) as 
there is only one non-trivial Lenard-Magri chain, we have to look for a single 
transversal vector field to deform the Poisson pencil Q — vP . 



16 



Let us consider 

Zo = Cl ^ + C "|? (5 - 1} 

and define 

X = QdH . (5.2) 

The vector field Z is the required ingredient for applying the bi-Hamiltonian 
setting to the Volterra lattice. This follows from the properties we list and 
prove below in a series of steps. 

a) Zq is a symmetry for the quadratic Poisson tensor P, 

Ue Zo P = 0. 

This follows noticing that the Jacobian of Z with respect to the coor- 
dinates (ci, . . . , Co) is, in terms of the standard generators Eij of n x n 
matrices, given by En + E nn , and by the explicit form of P. 

b) Still taking this property into account, one can easily verify that the 

action of Z on Q is given by 



where 



Lie Zo Q = Z A W , (5.3) 



9 9 \ d 9 

W = cic n ( — + cic 2 - c n _ic„- . (5.4) 

OCx OC n I OC2 Cc n _i 



Liez (Ho) = -H . (5.5) 



c) 

d) If n = 2k + 2, then 

Ue Zo TCk+i = 0, (5.6) 

where TCk+i is given by (|4.9|) . These last two properties can be easily 
verified by straightforward computations. 

Hence we can state 

Proposition 5.1 Let Z := — — — Z be the normalized symmetry of P. The 

Ho 

two bivectors 

P, Q = Q-Z AX (5.7) 



17 



form a Poisson pencil that restricts to the (generic) symplectic leaf S of P. 
The functions TCi are in involution also with respect to the deformed Poisson 
bracket {■,■}' associated with the bivector Q; hence their restrictions Tii define 
a bi-Lagrangian foliation of S. 

We are now left with the characterization of the Nijenhuis coordinates as- 
sociated with the restriction of the pencil Q — vP to the symplectic leaf 
S. 

Let us consider the vector field Z and the Hamiltonians Tii, with i — 1, . . . , k 
if n = 2k + 1, and % = 1, . . . , k + 1 if n = 2k + 2. 

Lemma 5.2 

Lie z (Hi) = Lie z ( -) = J % , (5.8) 

\7T J C!=c n =0 

that is, the Lie derivative ofTii with respect to Z is nothing but the numerator 
of Hi, evaluated at c± = c n = 0. 

Proof. The proof is a simple chain of computations. We report it here since 
this Lemma is crucial for the conclusion of the paper. 

We remark that, by the definition of the functions Jj and that of Z , we 
have 

Lie Zo (Ji) = Ji , (5.9) 

C2=-"=e n -i=0 



along with (eq (jS3J)) Lie^ (7Yo) = — ^o- Thus 

) = Ji 



Ue z (Hj) = - (^-Lie Zo (H Ji)) = -(Ji - J» 



=Cn_i=0 



Cl=Cn=0 

(5.10) 



□ 

From the Lemma above, it immediately follows 

Proposition 5.3 The polynomial Liez{H(fi)) factors as 

Lie z (H(\)) = v a A(v), (5.11) 

where a = 1 if n is odd, a = 2 if n is even, and A is a monic degree k 
polynomial in v = fi 2 which, thanks to A5.1U\) . is invariant along Z . 
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From the results recalled after Definition 13. 1\ we thus recover the eigenvalues 
of the Nijenhuis tensor on S as the roots of A(i/). Since they clearly are 
functionally independent, we can choose, as first half of Darboux-Nijenhuis 
coordinates, their square roots: 

Ci = V^- (5-12) 

Also, if we consider the normalized spectral curve relation, we see that the 
solutions Aj of the equation 

r(A,C0 = O (5.13) 

are invariant under Z as well, so that choosing one of the two solutions of this 
equations will provide a natural candidate for the remaining half of Nijenhuis 
coordinates. 

To show that actually this is the case, and simultaneously define a set of 
Darboux-Nijenhuis coordinates (that is, normalized Nijenhuis coordinates), 
we can rely, once more, on the results of |23*1 ITU] , thanks to the following 

Proposition 5.4 Let us consider the Lax operator £ of \4-<$ an d ^ ip be a 

Baker Akhiezer vector, for C, normalized with ipi — 1. Then on the k points 
Pi — (Aj, Ci) chosen according to the above recipe, if) has a pole. 

Proof. The normalized BA function ip has poles in the zeros of the (1,1) 
element of the classical adjoint matrix 

{H-Cy. (5.14) 

This is the determinant of the (n — 1) x {n — 1) matrix 

n-2 

M = 111 - ^ ^C a +l(E a ,a+l + E a +l,a), (5.15) 
a=l 

whose determinant equals the A(^z 2 ) if n is odd, and fxA(fx 2 ) if n is even. 

□ 

This shows that the functions (Ci, \)i=i,...,k selected according the bi-Ham- 
iltonian scheme herewith presented do indeed coincide with those found by 
Veselov and Penskoi via the method of poles of the BA function. 

As a closing remark, we notice that a set of canonical coordinates for 
the Volterra lattice can be obtained via "purely bi-Hamiltonian methods" as 
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follows. As it has been remarked in Proposition 12. 2\ a possible path is to 
use the Hamiltonian vector fields associated with the coefficients Pi of the 
minimal polynomial of the Nijenhuis tensor A (A) to deform the polynomial 
Casimir of the Poisson pencil. In this way we obtain new polynomials that 
satisfy the characteristic equation of a Nijenhuis functions generator ([2.4)1 . 
and hence we can use them to generate Darboux-Nijenhuis coordinates. 

For the sake of concreteness we will stick to the case of an odd number 
of sites n. 

We consider the minimal polynomial of the induced Nijenhuis tensor, 
given by (according to Proposition 15.31 with n — 2k + 1), 



A(i/) = (-Ue z (H(fi))) = v k -p x v k 1 p k , where pj = -lAe z (Hj). 

V \i J fl 2 =U 

(5.16) 

Thanks to the explicit characterizations of the Hamiltonians Tii and of the 
vector field Z, and taking Lemma 15.21 into account, it is not difficult to 
ascertain that 

k 

Pk = Y[c 2i . (5.17) 

i=l 

Keeping into account the explicit form of the quadratic Poisson tensor ()4.3|) . 
the Hamiltonian vector field associated with \ogpk is given by the very simple 
expression 

Y = Pd\og Pk = Cl A_ Cj A (5.i 8 ) 
act oc n 

If we define as "first half of the Darboux-Nijenhuis coordinates the loga- 
rithms of the eigenvalues of the Nijenhuis tensor, rather than their square 
roots, that is, if we consider 

(pi = l0g(l/;), 

we clearly have that, in terms of the canonical coordinates ipi conjugated to 
the <f>i we are seeking, it holds 

Y = Pj2d (j)l = J2^-, (5-19) 

i i 

so that according to the recipe we are using, we need to find an exact eigen- 
function generator satisfying 

Y(V(u)) Elmod A(z/). 
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Thanks to the explicit form of the vector field Y we can easily establish, 
arguing as in the proof of Lemma 15.21 the following equalities: 

J'i ■= Liey(Ji) = Jj(ci, 0, . . . , 0, -c n ), 

J'! := Liey(j;) = J<(ci, 0, . . . , 0, c n ), (5.20) 

and, finally, J-" = J[. 

This shows that the functions 

rn = log(^), A = log (y(J2 ^7i k ) + Y\Y1 U * H ^ ) , * = 1, • • • , fc- 

\ fc k J 

provide a set of canonical DN coordinates for the Volterra lattice with 2k + 1 
sites. An analogous result holds for the VL with an even number of sites, 
although the vector field Y has a more complicated expression. 
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